UNCLASSIFIED 


_ AD  NUMBER _ 

AD841124 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
and  their  contractors;  Critical  Technology;  SEP 
1968.  Other  requests  shall  be  referred  to  Rome 
Air  Development  Center,  Attn  EMASA,  Griffins 
AFB,  NY  13550.  This  document  contains  export- 
controlled  technical  data. 


_ AUTHORITY 

RADC  ltr,  14  Oct  1971 


THIS  PAGE  IS  UNCLASSIFIED 


■ 


When  US  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose  other 
than  a  definitely  related  government  procurement  operation,  the  government  thereby  incurs 
no  responsibility  nor  any  obligation  whatsoever;  and  the  bet  that  the  government  may  have 
formulated,  furnished,  or  in  any  way  supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded,  by  implication  or  otherwise,  as  in  any  manner  licensing  the 
holder  or  any  othe»  person  or  corporation,  or  conveying  any  rights  or  permission  to  manu¬ 
facturer,  use,  or  sell  acy  patented  invention  that  may  in  any  way  be  related  thereto. 


accessiwj* 


IUHWN3«(IC£B 

JUSTtFlCtttW 


niTitiBinioa/avaitABiuiTt  ca 

BiTT.  I  AVAIL  and/or  SfKUt 


Do  not  return  this  copy.  Retain  or  destroy 


Best 

Available 

Copy 


*  immmi  iW* 


HF  ANTENNA  GAIN  PATTERN  MEASUREMENTS: 
EXPANDED  LITTLE  IDA  PROGRAM 

Joseph  J.  Simons 
Frederick  C.  Wilson 
Richard  A.  Schneible 


This  document  is  subject  to  special 
export  controls  and  each  transmittal 
to  foreign  governments,  foreign  na¬ 
tionals  or  representatives  thereto 
may  be  made  only  with  prior  approval 
of  RADC  (EM ASA),  GAFB,  N.Y.  13440. 


FOREWORD 


This  report  was  prepared  under  Project  5582,  Task  558202,  by  personnel  of 
Rome  Air  Development  Center,  Environmental  Studies  Section,  EMASA,  The  RADC 
Flight  Test  Division  (EMF),  the  Measurement  Section  and  Test  Operations  Unit 
(EMCVM)  of  the  RADC  Communications  Division,  and  the  Heavy  Military  Electronics 
Division  of  General  Electric  Company,  Syracuse,  N.  Y,  ,  assisted  in  obtaining  the 
data  discussed  herein.  Specif!  acknowledgment  is  given  to  Donald  McCoy  and  Robert 
Wengenroth  of  General  Electric;  and  to  Donald  Long,  Charles  Lupica,  Capt.  Robert 
Welch,  Capt.  William  Freeman,  and  Capt.  James  Hayes  of  RADC,  EMCVM. 

Distribution  of  this  report  is  restricted  under  the  U.S.  Mutual  Security  Acts  of 
1949  to  protect  techniques  revealed  herein. 

This  Technical  Report  has  been  reviewed  and  is  approved. 


Approved: 


Acting  Chief,  Space  Surveillance 
and  Instrumentation  Branch 
Surveillance  5  Control  Division 


t 


ii 


ABSTR  k  " 


This  report  describes  the  data  collection,  reduction,  and  processing  techniques 
employed  in  an  antenna  gain  measurements  program.  This  gain  measurements  pro¬ 
gram  is  an  integral  part  of  the  overall  Expanded  Little  IDA  effort. 

The  gain  patterns  of  transmitting  antennas  located  at  Thule,  Greenland; 
Keflavik,  Iceland;  and  Coco  Solo,  Panama  Canal  Zone  were  determined  and  the  re¬ 
sults  reported. 

Site  and  aircraft  instrumentation  used  to  obtain  the  data  and  manual  and  com¬ 
puter  processing  methods  used  to  reduce  the  data  are  discussed. 
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I.  INTRODUCTION 


This  report  describes  the  measurement  and  data  processing  techniques  used  in 
determining  the  gain  patterns  of  full-scale  high  frequency  (HF)  antennas  located  at 
four  sites  utilized  in  the  RADC  Expanded  Little  IDA  Program  (ELIDA).  Transmitting 
sites  are  located  at  Pingorssuit  "P"  Mountain,  Thule  AFB,  Greenland;  Site  Dye  5, 
Keflavik  NAS,  Iceland;  and  Coco  Solo  Test  Annex,  Panama  Canal  Zone.  The  transmit 
sites  beam  their  signals  to  a  central  receiving  site  located  at  the  Starr  Hill  Test  Annex, 
Remsen,  N.Y. ,  situated  19  miles  north  of  Griffiss  Air  Force  Base. 

The  measurement  technique  involved  the  utilization  of  an  Air  Force  NKC-135  re¬ 
fueling  tanker  instrumented  with  several  KF  receiving  chains.  The  data  processing 
and  reduction  procedures  involved  both  manual  and  computer  techniques. 

The  antenna  measurements  were  performed  in  support  of  contract  AF30(602)- 
3946,  Expanded  Little  IDA,  with  the  General  Electric  Company,  Heany  Military  Elec¬ 
tronics  Division,  Syracuse,  New  York.  GE  operated  and  maintained  these  sites  under 
this  effort.  Both  RADC  and  GE  personnel  participated  in  the  measurement  program. 
Data  reduction  was  primarily  performed  by  the  Environmental  Studies  Section  of 
RADC,  with  the  consultant  services  of  GE  engineers. 

This  report  discusses  the  techniques  employed  and  results  obtained  for  the 
antenna  gain/flight  tests  measurements  of  the  antennas  located  at  the  Thule,  Iceland, 
and  Panama  transmitting  sites.  The  Panama  and  Thule  tests  were  conducted  during 
August  and  October  1966,  while  the  Iceland  measurements  were  conducted  during 
October  1967.  Some  preliminary  measurements  have  been  made  at  the  Starr  Hill  re¬ 
ceiving  site,  but  additional  tests  are  required.  Results  of  Starr  Hill  measurements 
will  be  reported  on  at  t  later  date. 

A  discussion  of  the  Expanded  Little  IDA  Program,  which  generated  the  re¬ 
quirement  for  antenna  gain  patterns,  is  presented  in  Appendix  A.  One  investiga¬ 
tion  involves  the  measurement  of  high  frequency  signal  losses  associated  with  the 
individual  modes  that  propagate  over  each  path.  The  analysis  of  propagation  loss 
requires  (a)  a  determination  of  the  total  loss  between  transmitter  and  receiver, 
and  (b)  the  separation  of  known  losses  and  gains. 

The  transmitting  and  receiving  antenna  gains  must  be  removed  from  the  loss 
measurements  so  that  the  data  will  be  independent  of  the  antenna  used.  This  ex¬ 
periment  is  being  performed  over  paths  that  have  a  2000  nm.  range,  which  is 
the  limiting  distance  for  propagation  via  the  one-hop  F  (1F2)  mode.  When  prop¬ 
agating  over  these  -anges  the  1F2  mode  is  received  at  elevation  angles  between 
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-1°  to  4°,  while  the  two-hop  F  mode  is  received  at  elevation  angles  between  12°  and 
20°.  Therefore,  to  perform  a  loss  measurement  on  the  individual  modes  of  propaga¬ 
tion,  the  antenna  gain  patterns  must  be  measured  as  a  function  of  elevation  angles  in 
the  range  from  -1°  to  +20°. 
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II.  SITE  INSTRUMENTATION 


A.  EQUIPMENT 

The  measurement  technique  employed  in  the  tests  of  the  Thule,  Iceland,  and 
Panama  transmitting  antennas  was  identical.  Signals  were  transmitted  from  the 
ground  and  received  in  the  NKC-135  aircraft.  A  block  diagram  of  the  inrv  unr.entation 
used  at  the  transmitting  sites  for  the  tests  is  shown  in  Figure  1. 

Four  transmitting  systems,  each  tuned  to  a  different  frequency,  weiv  time 
multiplexed  to  the  antenna  under  test.  Each  transmitting  system  consisted  of  (1)  a 
frequency  synthesizer,  (2)  an  RF  exciter,  and  (3)  a  1  kw  to  3  kw  amplifier.  These 
systems  were  sequentially  connected  in  turn  to  the  antenna  under  test  by  a  sequence 
keyer.  Each  transmitter  was  keyed  on  for  one  second  and  off  for  three  seconds. 
Therefore,  four  frequencies  were  transmitted  in  four  seconds.  Power  monitoring 
was  accomplished  by  connecting  an  oscilloscope  to  a  directional  coupler  so  as  to  main¬ 
tain  a  constant  power  level  for  all  transmitters  during  each  of  the  tests. 

The  pulsing  rate  of  one  second  on,  three  seconds  off,  provided  adequate  sam¬ 
pling  of  the  antenna  pattern,  permitted  four -frequency  operation  without  complex 
multicoupler  techniques,  and  allowed  the  noise  and  interference  levels  to  be  checked 
between  pattern  samples.  Thus,  in  effect,  four  gain  patterns  were  measured  in  one 
radial  aircraft  flight  pattern,  reducing  flying  time  by  75  percent.  The  data  rate  using 
this  pulsing  technique  was  adequate  for  measuring  patterns  below  about  25°  elevation 
angle.  Above  20°  the  data  rate  is  too  slow  to  represent  the  fast -changing  pattern. 

The  ability  to  measure  the  noise  and  interference  level  permitted  the  analysis  of 
some  data  which  would  otherwise  have  been  unusable  because  it  indicated  when  noise 
levels  were  high  and  also  aided  in  selecting  frequencies  that  correspond  to  low  noise 
levels.  Also,  when  signal  levels  are  close  to  the  noise  level,  a  separate  measure¬ 
ment  makes  it  possible  to  correct  for  its  influence  on  the  signal. 

B.  ANTENNAS 
1.  Thule 

A  Granger  Associates,  Model  726-5,  antenna  was  installed  at  the  MP"  Mountain 
Site.  This  antenna  is  located  or  a  rocky  point  between  two  troposcatter  antenna  in¬ 
stallations.  From  this  site  there  is  an  unobstructed  path  south,  with  the  horizon  be¬ 
ing  Baffin  Bay  2600  feet  below.  The  dip  angle  to  the  horizon  is  about  1°  from  the 
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Figure  1.  Transmitting  Equipment 


2600  ft.  elevation  of  the  site.  Figures  2  and  3  are  photos  of  the  Thule  antenna.  The 
specifications  of  the  Granger  726-5  are: 

Frequency  Range  -  4-64  MHz 

Type  -  Log -Periodic  Mono  pole 

Application  -  Receive  and  Transmit 

Power  Capacity  -  1  KW  average,  30  KW  PEP 

Input  Impedance  -  50  ohms 

Polarization  -  Vertical 

2.  Keflavik  (Dye  5),  Iceland 

A  Granger  726-5  LPV  (identical  to  that  at  "PM  Mt. )  was  also  installed  at  Site 
Dye  5,  Iceland.  The  Iceland  antenna  is  located  at  a  seaside  site,  very  close  to  the 
ocean.  The  patterns  of  this  antenna  were  subject  to  tidal  effects. 

3.  Coco  Solo,  Panama  Canal  Zone 

The  measurement  antenna  at  Coco  Solo  is  a  log -periodic,  vertically  polarized 
dipole,  ALL-Products  Co.  Model  LPV  12-A.  This  antenna  also  has  a  frequency  range 
of  4  to  64  MHz  and  its  specifications  are  similar  to  the  Thule  antenna.  The  Coco  Solo 
site,  located  on  the  Atlantic  side  of  the  Canal,  is  situated  in  a  swampy  area  about  one 
foot  above  sea  level.  The  ocean  tides  have  a  range  of  only  about  nine  inches.  The 
data  have  shown  that  this  small  tide  has  no  measurable  influence  on  the  antenna  pattern. 
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Figure  2.  Thule  Antenna 


Figure  3.  Thule  Site  -  Looking  Along  Bo  re  sight 


III.  AIRCRAFT  INSTRUMENTATION 


An  NKC-135  RADC  aircraft  was  instrumented  with  (1)  a  four-channel  receiving 
and  recording  system  which  is  compatible  with  the  high  frequency  transmitting  system 
and  (2)  a  receiving  loop  antenna  which  could  be  oriented  to  receive  either  vertically 
polarized  or  horizontally  polarized  radiation.  A  block  diagram  of  the  aircraft  instru¬ 
mentation  is  shown  in  Figure  4.  Photos  of  the  aircraft  and  receiving  loop  are  shown 
in  Figures  5  and  6. 


Once  each  signal  has  been  received  by  the  two -turn  loop  *ud  balun,  it  is  ampli¬ 
fied  and  passed  through  a  switching  attenuator/directional  coupler  network.  As  illus¬ 
trated  on  the  diagram,  the  signal  is  then  coupled  to  one  of  four  receiving  and  recording 
chains.  Each  chain  is  composed  of  a  step  attenuator,  and  an  R-390  HF  receiver.  The 
AGC  voltage  level  from  each  receiver,  which  is  nearly  logarithmic  with  signal  level, 
is  recorded  on  a  four  channel  Brush/Mark -200  analog  recorder.  As  shown  in  Figure 
4,  a  Hewlett-Packard  Frequency  Synthesizer,  Model  No.  5100/5110,  and  an  attenuator 
are  used  in  calibrating  each  receiver  channel. 


8 


Aircraft  Receiving  Instrumentation 


Figure  5.  NKC-135  Flight  Test  Aircraft 
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Figure  6.  Receiving  Loop 


IV.  DATA  COLLECTION  PROCEDURE 


The  same  technique  was  used  at  all  locations  for  the  collection  of  antenna  pat¬ 
tern  data.  At  each  site  and  for  each  antenna  gain  measurement,  the  aircraft  flew  a 
radial  flight  plan  out  along  the  boresight  of  the  antenna  and  also  at  plus  and  minus  5° 
and  10°  off -bo resight.  Data  taking  commenced  when  the  aircraft  was  directly  over  the 
site  (at  the  zenith  of  the  antenna).  The  radial  patterns  were  flown  to  ranges  of  300 
nautical  miles  from  the  site,  and  at  altitudes  between  27,000  and  31, 000  feet.  How¬ 
ever,  altitude  was  kept  constant  through  each  individual  run.  With  respect  to  antenna 
elevation  angles,  these  radial  flights  correspond  to  elevation  angles  of  90°  to  -1°.  As 
previously  mentioned,  primary  interest  was  centered  on  data  taken  between  elevation 
angle?  of  -1°  and  20°,  which  covers  the  range  over  which  the  lower  order  dominant 
modes  (IF,  2E,  2F,  etc.)  are  launched. 

The  gain  pattern  of  each  antenna  was  measured  as  a  function  of  eight  frequencies 
in  the  9  to  30  MHz  region.  Therefore,  a  total  of  forty  (8  frequencies,  5  radials  per 
frequency)  successful  data -taking  radials  were  required.  To  check  the  validity  of  the 
data,  the  boresight  radials  were  repeated.  In  addition,  radial  missions  were  sched¬ 
uled  to  check  the  cross-polarized  response  of  the  antenna  under  test.  This  was  ac¬ 
complished  by  changing  the  polarization  of  the  aircraft  loop  antenna  (i.  e. ,  from 
vertical  to  horizontal  polarization). 

Range  information  (that  is,  knowing  the  location  of  the  aircraft  at  any  given 
time)  is  one  source  of  error  in  gain  measurements.  Navigation  procedures  were  set 
up  so  that  the  range  of  the  aircraft  from  the  transmitting  antenna  could  be  determined 
continuously  as  the  aircraft  flew  along  a  specific  heading.  For  an  aircraft  flying  at  a 
constant  altitude,  the  location  or  specific  range  of  the  aircraft  from  the  antenna  under 
test  corresponds  to  a  specific  antenna  elevation  angle. 

In  the  case  of  the  Panama  and  Iceland  measurements,  a  doppler -radial  navigation 
radar,  APN-81,  located  on  the  aircraft  was  used  in  determining  range.  Range  marks, 
in  five  and  ten  nautical  mile  increments,  were  recorded  on  the  same  recorders  used 
for  the  signal  test  data.  This  navigational  syslem  is  specified  as  having  a  range  ac¬ 
curacy  of  plus  or  minus  one  nautical  mile.  Fo.r  the  Thule  measurements,  this  system 
was  supplemented  by  use  of  the  approach  radar,  FPS-47,  located  at  the  main  base  at 
Thule.  This  radar  was  used  to  position  the  aircraft  directly  over  the  site  on  each 
radial  run.  In  tho  Panama  and  Iceland  tests,  visual  techniques  were  employed  in  posi¬ 
tioning  the  aircraft  so  as  to  establish  the  site  as  zero  reference  for  each  radial. 
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V.  DATA  REDUCTION 


A.  THE  LOSS  EQUATION 

The  signal  power,  pr,  present  at  the  receive  antenna  terminals  of  the  flight  test 
aircraft  is  given  by  the  familiar  equation: 


P 

r 


Wr 


where: 


p^  -  power  output  at  the  terminals  of  the  transmitter 

g^  =  gain  of  the  transmitting  antenna 

g  =  gain  of  the  receiving  antenna 
r 

1,  _  =  free  space  loss 

bi 

equation  (1)  can  be  expressed  in  decibels  as 


P 

r 


gt  + 


(1) 


(2) 


By  rearranging  terms,  the  gain  of  the  transmitting  antenna  can  be  expressed  as: 

gt=Pr-Vgr+  Srf  (3> 

Thus  the  gain  of  the  transmitting  antenna  can  be  calculated  provided  the  transmitted 
power,  received  power,  receiving  antenna  gain,  and  free  space  loss  are  known. 

1.  Transmitted  Power 

Power  levels  were  determined  at  the  transmitter  site  by  using  a  Tektronix  oscil¬ 
loscope  to  measure  the  voltage  across  a  CO-ohm  load  in  the  arm  of  a  Sierra  directional 
coupler  connected  between  tl'2  transmitter  and  the  antenna.  At  the  end  of  the  program, 
the  oscilloscope  was  checked  by  inserting  a  20  dbm  signal  from  an  HP  606A  signal 
generator.  Power  was  monitored  throughout  each  data  run  so  that  any  lev  >1  changes 
could  be  noted. 
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Received  Power 


On  the  aircraft,  received  power  is  recorded  directly  on  chart  recorders  using 
calibrated  receivers.  Calibration  run  are  made  for  each  radial  to  check  drift.  At¬ 
tenuators  are  inserted  in  the  receiver  chain  so  that  the  recorders  are  kept  within  their 
linear  range  throughout  the  data  period.  The  value  of  the  attenuators  in  the  circuit 
along  with  range  marks  from  the  navigator  are  entered  on  the  chart  records  by  the 
operator. 

3.  Loop  Gain  -  Gr 

The  calculation  of  the  loop  gain  is  one  of  the  most  difficult  parts  of  this  experi¬ 
ment.  In  order  to  calculate  this  correctly  it  is  necessary  to  fly  the  aircraft  against 
an  antenna  whose  pattern  is  known.  A  horizontal  dipole,  located  at  RADC's  Verona 
Test  Site  will  be  used  to  calibrate  the  loop  for  horizontal  polarization.  Measurement 
techniques  to  determine  the  gain  of  the  loop  for  vertical  polarization  are  being  investi¬ 
gated.  This  investigation  involves  a  determination  of  a  standard  reference  antenna 
whose  pattern  is  well  known.  An  interim  method  was  devised  that  involvod  the  selec¬ 
tion  of  a  point  on  the  antenna  pattern  at  an  elevation  angle  that  was  relatively  free  of 
ground  effects.  The  gain  at  this  point  was  set  equal  to  the  theoretical  gain  of  the  an¬ 
tenna.  Since  the  relative  gains  were  known,  this  gave  the  absolute  antenna  patterns. 


4.  Free  Space  Loss 


The  free  space  loss  incurred  by  a  high  frequency  wave  of  wavelength  X  trans¬ 
fer  sing  a  distance  R  is  given  fay: 


(4) 


In  our  case,  R  is  the  slant  range  from  the  transmitting  site  to  the  airplane  and  is 
calculated  in  the  following  maimer.  Referring  to  Figure  7 

R 

0  (radians)  =—^~  (5) 

and  fay  the  law  of  cosines 

R2  =  a2  +  (a+-h)2  -  2a(a+h)  cos  0  (6) 

Since  X  is  known,  free  space  loss  can  be  calculated  knowing  the  height  of  the  airplane 
and  the  ground  range  between  the  transmitter  and  airplane. 


TRANSMITTER  POSITION 
AIRPLANE  POSITION 
RADIUS  OF  EARTH 
HEIGHT  OF  AIRPLANE 
SLANT  RANGE  TO  AIRPLANE 
GROUND  RANGE  TO  AIRPLANE 
EARTH’S  CENTRAL  ANGLE 
TAKEOFF  ANGLE 


Figure  7.  Test  Geometry 
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B.  ELEVATION  ANGLE 

1.  Geometrical  Angle 

Although  elevation  angle,  as  such,  does  not  appear  in  the  basic  equation  for  de¬ 
termining  gain  of  the  transmitting  antenna,  it  is  required  for  data  presentation  so  that 
the  antenna  gain  patterns  can  be  made  as  a  function  of  a  common  factor.  With  gains 
plotted  versus  elevation  angle,  direct  comparisons  of  the  gain  change  with  frequency, 
off-bo  re  sight  angle,  etc. ,  can  be  made.  This  method  of  data  presentation  is  advan¬ 
tageous  in  comparing  the  patterns  of  different  antennas. 

Referring  again  to  Figure  7,  the  geometrical  angle  &q  can  be  calculated  by  using 
the  law  of  cosines: 


2  2  2 

(a+h)  =  R  +  a  -  2aR  cos  a, 

(7) 

or 

cos  a  =  (R2  +  a2  -  (a*h)2)/ : 

(8) 

and  since  a  =  90  + 

(9) 

cos  (90  +  <*0)  =  T*2  +  a2  -  (a+h)2/2aR, 

(10) 

and 

-sin  aQ  =  (R2  +  a2  -  (a2  +  2 ah  +  h2))/2aR, 

(11) 

or 

sin«0  =  (2 ah  +  h2  -  R2)/2a R 

(12) 

All  the  quantities  in  the  above  equation  are  known  and  thus  the  calculation  of  the 
geometrical  elevation  angle  can  be  made. 

2.  Tropospheric  Correction 

Due  to  tropospheric  refraction,  the  geometrical  angle  is  not  the  angle  at  which 
the  energy  is  radiated  from  the  antenna.  It  is  necessary  to  calculate  the  antenna  pat¬ 
terns  in  terms  of  this  latter  angle  which  we  shall  call  the  true  angle.  This  can  be  seen 
from  Figure  8,  which  shows  how  the  energy  from  a  transmitter  located  at  T,  leaves 
the  antenna  at  different  angles  to  intercept  two  receivers,  Ri  and  R2,  located  at  the 
same  geometrical  angle.  Thus  antenna  patterns  flown  at  different  heights  would  give 
different  results  unless  these  refinements  were  made.  These  corrections  in  elevation 
angle  must  be  made  before  the  data  is  presented,  and  have  been  made  part  of  the 
computer  program,  ANTPAT,  which  is  discussed  in  the  next  section.  The  computer 
employs  a  table  look-up  technique  using  values  for  the  refraction  data  as  published  in 
Reference  7. 
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T  «  TRANSMITTER  SITE 
Rl  =  RECEIVER  #1 
R9  -  RECEIVER  #2 

§  -  GEOMETRICAL  ANGLE  TO  R-j  and  R? 
0,  -  TRUE  ANGLE  TO  Ri 
02  -  TRUE  ANGLE  TO  R2 


Figure  8.  Refraction  Effects 
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VI.  COMPUTER  PROCESSING 


For  efficient  processing  of  the  antenna  pattern  data,  three  basic  computer  pro¬ 
grams  wer*»  written  in  Fortran  for  use  with  R ADC's  Control  Data  1604  Computer. 
Programs  written  were  named  ANTPAT,  TAPPLT,  and  DATAVE.  Listings  of  these 
programs  can  be  found  in  Appendix  B.  Comment  cards  have  been  added  to  identify  the 
main  sections.  Sample  outputs  of  each  of  the  programs  can  be  found  in  Section  VII  of 
this  report. 

A.  PROGRAM  ANTPAT 

Figure  9  is  a  block  diagram  of  Program  ANTPAT.  ANTPAT  is  the  basic  data 
reduction  program.  It  uses  the  antenna  pattern  data  which  is  punched  onto  input  cards, 
performs  all  the  necessary  calculations,  i.  e.  ,  slant  range,  free  space  loss,  etc. ,  and 
provides  a  tabulation  of  antenna  gain,  take-off  angle,  and  other  selected  variables.  It 
records  all  information  necessary  for  further  processing  by  the  other  programs  on  a 
magnetic  tape. 

1.  Data  Requirements 

The  information  necessary  for  ANTPAT  to  calculate  the  gain  of  the  transmitting 
antenna  can  be  summarized  from  Section  V  as  follows: 

•  transmitter  power 

•  received  power 

•  attenuator  settings 

•  receiver  calibration  data 

•  frequency  of  transmission 

•  height  of  aircraft 

•  ground  range  from  transmitter  to  aircraft 

•  loop  gain 

•  tropospheric  angle  corrections 

In  addition  certain  information  is  necessary  to  identify  features  of  the  data  run 

itself: 

•  site  identity  and  loop  polarization 

•  date  of  the  run 
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CALCULATE  RECEIVED  SIGNAL  STRENGTH 
USING  THE  CALIBRATION  CURVE 


Figure  9.  Block  Diagram  -  Program  ANTPAT 
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MW*1 


•  heading 

•  run  number  (the  sequential  number  of  the  run  in  the  flight  test  plan) 

•  time  of  day  (needed  to  help  assess  certain  geophysical  effects,  e. g. , 
weather,  tides,  etc.) 

•  identifying  number  of  the  receiver  used  for  a  specific  run  and  frequency 

The  tropospheric  corrections  and  loop  gain  calculations  were  made  an  integral 
part  of  ANT  PAT  and  are  not  a  part  of  the  input  data.  Transmitter  power  was  meas¬ 
ured  at  the  transmit  site.  Received  powers  and  calibration  data  were  recorded  by  the 
receiving  system.  From  a  practical  point  of  view,  it  was  found  to  be  a  difficult  and 
time-consuming  task  to  convert  the  data  points  directly  to  power  received  in  dB*a 
using  the  calibration  curve.  For  each  run,  therefore,  the  points  on  the  calibration 
curve  and  the  data  points  were  read  with  reference  to  the  divisions  on  the  chart  paper. 
The  computer  calculates  the  data  point  locations  on  the  calibration  curve  and  their 
corresr mding  received  power  values. 

Readings  of  the  data  points  were  mrde  usually  in  five  nautical  mile  intervals  out 
to  60  miles  and  every  ten  nautical  miles  thereafter.  Any  obvious  peaks  and  nulls  in 
the  data  which  do  not  occur  at  these  specific  range  intervals  were  recorded  so  that 
their  full  extent  could  be  assessed. 

The  site  identity  and  polarization  are  defined  in  the  computer  program  by  a  vari¬ 
able  name,  NPOL.  This  variable  must  be  specified  by  a  Fortran  statement  before 
each  computer  run.  Consequently  every  computer  run  must  consist  of  data  from  the 
same  site  and  of  the  same  polarization. 

All  of  the  other  quantities  listed  above  are  either  entered  as  part  of  the  Flight 
Test  log  or  as  notes  on  the  chart  recorder  paper. 

2.  Input  Data  Format 

ANTPAT  was  written  as  an  input  program  for  all  the  data  specified  in  Section 
VL  A.1  above.  Data  must  be  supplied  to  the  program  in  the  following  order  and 
format. 

a.  Information  header  card 


Columns 

Format  * 

Variable 

1-6 

16 

date 

7-9 

13 

heading 

10-11 

12 

run  number 

21 

Columns 

Format 

Variable 

12-15 

14 

time 

16-25 

FI  0.0 

altitude  (K  feet) 

Data  header  card 

Columns 

Format 

Variable 

1-2 

12 

number  of  data  points 

3-4 

re 

receiver  number 

5-14 

F10.0 

frequency  (MHz) 

Power  card 

Columns 

Format 

Variable 

1-10 

F10.0 

power  in  DBM 

Calibration  Data 

Columns 

Format 

Variable 

1-4,  5-8, 

— ,69-72 

18F4.1 

calibration  data  (DIV) 

Data  Points 

Columns 

Format 

Variable 

1-5,  6-11, 
-,71-75 

15F5. 2 

ground  range  (NM) 

1-5,  6-11, 

— ,71-75 

15F5. 2 

data  points  (DIV) 

1-5,  6-11, 
-,71-75 

15F5. 2 

attenuations  (DB) 

3.  Calculations  and  Output  Form 

Attar  reading  the  required  input  data  as  above,  ANTPAT  makes  all  essential 
calculations  as  described  in  Section  V.  The  determination  of  loop  gain  is  done  sep¬ 
arately  in  a  subroutine.  The  other  calculations  are  all  straight-forward  and  will  not 
be  elaborated  on  further.  The  block  diagram  (Figure  9)  and  the  listing  in  Appendix  B 
should  enable  the  reader  to  follow  the  program  with  little  difficulty. 


The  tabular  output  of  ANTPAT  (see  Section  VII  for  an  example)  provides  a  quick 
look  at  the  data  to  check  for  obvious  errors  before  further  data  processing  is  carried 
out.  In  addition,  ANTPAT  writes  all  data  necessary  for  the  subsequent  programs  on 
a  magnetic  tape.  All  data  points  are  included  in  the  tabular  output  of  ANTPAT.  Data 
taken  at  angles  greater  than  twenty-five  degrees  are  of  a  secondary  interest  and  are 
not  written  on  the  tape.  The  first  time  the  program  is  run,  Section  A  (see  Figure  9) 
is  removed.  Data  are  recorded  on  tape  followed  by  a  dummy  record.  For  each  addi¬ 
tional  run,  Section  A  is  reinserted.  The  program  adds  on  to  the  data  tape,  erasing 
the  dummy  record  and  moving  it  to  the  end.  The  dummy  record  serves  as  a  flag  for 
subsequent  programs  to  indicate  the  end  of  the  data, 

B.  PROGRAM  TAPPLT 

Program  TAPPLT  is  the  basic  data  plotting  program.  It  uses  the  tape  written 
by  ANTPAT  and  performs  a  search  for  selected  data.  It  then  writes  a  magnetic  tape 
which  is  compatible  with  RADC's  EAI  3500  data  plotter  system.  Gain  versus  elevalioD 
angle  plots  can  be  obtained  off  line  using  this  system.  Figure  10  ir>  a  basi.  block  dia¬ 
gram  of  TAPPLT. 

TAPPLT  plots  the  repeat  runs  for  a  given  site,  frequency,  heading,  and  polar¬ 
ization  on  the  same  graph  so  that  the  repeatability  of  the  patterns  as  wrell  as  the  pat¬ 
terns  themselves  can  be  studied.  A  major  part  of  the  program  deals  with  the  output 
of  header  information  and  the  drawing  and  labeling  of  the  horizontal  and  vertical  axes. 
Many  statements  are  calls  to  the  CDC  1604  machine  subroutines  and  the  programming 
of  these  statements  will  not  be  discussed  here.  A  detailed  discussion  of  these  sub¬ 
routines  can  be  found  in  Reference  9. 

C.  PROGRAM  DATAVE 

Program  DATAVE  also  utilizes  the  tape  generated  by  ANTPAT.  It  calculates 
average  gain  and  standard  deviation  of  the  data  points.  Similar  to  TAPPLT,  it  writes 
a  tape  for  use  with  the  EAI  data  plotter. 

Figure  11  is  a  basic  block  diagram  of  DATAVE.  It  consists  of  three  separate 
parts:  the  main  program,  subroutine  AVERAGE,  and  subroutine  PLOT. 

The  main  program  searches  the  tape  for  selected  data.  AVERAGE  then  takes 
repeat  runs  for  a  given  site,  frequency,  heading,  and  polarization  and  examines  them 
at  common  values  of  elevation  angles.  At  these  points  it  calculates  their  average 
gain  and  standard  deviation. 

1.  Average  and  Standard  Deviation 

The  calculation  of  average  and  standard  deviation  requires  that  we  must  work 
with  gain  itself  and  not  the  logarithm  of  gain.  The  statistical  analysis  can  then  be 
performed.  Afterwards  the  average  and  standard  deviations  can  be  converted  back 
to  decibels. 
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START 


PRELIMINARY  CALCULATIONS 


SET  PARAMETERS 
FOR  SEARCH 


SEARCH  COMPLETE? 


REWIND  TAPE 
NGO  -  1 


Figure  10.  Block  Diagram  -  Program  TAPPLT 
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Figure  11.  Block  Diagram  -  Program  DATAVE  Main  Program 
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COMMON  DATA  DIMENSION 


PUT  CURVES  IN  ORDER 
(BY  ANGLES) 


FIND  THE  LONGEST  CURVE 


CALCULATE  GAIN  FOR  EACH 
CURVE  AT  THIS  POINT 


SET  STANDARD  DEVIATIONS 
EQUAL  TO  ZERO 


REPEATED 
FOR  EACH 
POINT 


CALCULATE  AVERAGE  GAIN 


CALCULATE  STANDARD  DEVIATION 


PRINT  AVERAGE  GAIN,  ANGLE, 
STANDARD  DEVIATIONS 


ONLY 
I  CURVE 


CALL  PLOT 


RETURN 


Figure  11.  Block  Diagram  -  Program  DATAVE  Subroutine  AVERAGE 
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AVERAGE  converts  the  gain  data  from  dB’s  to  multiplicative  factors  (e.  g.  ,  3 
dB  is  converted  to  2).  The  program  then  calculates  the  average  (X)  and  standard 
deviation  (S*)  of  these  multiplicative  factors.  The  average  in  db  (?)  is  calculated 
using  the  relationship: 


Y  =  10  log  X  •  (13) 

The  formula  for  the  standard  deviation  in  dB  (Sy)  can  be  derived  from  the  definition: 


thus 


N  2  N  2 
5Z  <y i-y)  =  E  0yt) 
s  2  1=1  i^l _ L 

y  N  N 


I  (» ‘t 


s2  = 
y 


N 

N 


2  £  <*Xi> 

.  2  /dY\2  i=i 

Sy  =\^x) - N - 


(14) 


(15) 


(16) 


but  for  the  dB  relationship 


y  =  10  log  X 


(17) 


dy  10 
dx  X 


therefore  we  obtain 


2.  Output  Form 


(18) 


(lf») 


Subroutine  AVERAGE  provides  a  tabulation  of  the  average  gain  and  standard 
deviation  versus  elevation  angle  (see  Section  VII  for  an  example). 

Subroutine  PLOT  writes  a  tape  for  use  with  the  Data  Plotter  system  similar  to 
TAPPLT  to  obtain  plots  of  average  gain  versus  elevation  angle.  All  headings  for  a 
specific  frequency  and  polarization  are  plotted  on  the  same  graph  so  that  the  gain  as 
a  function  of  angle  off-boresight  can  be  studied  easily.  An  example  of  this  plot  can  be 
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found  in  Section  VTL  The  standard  deviation  cannot  be  calculated  for  azimuths  where 
data  was  recorded  only  once.  In  that  case,  the  standard  deviation  is  set  equal  to  zero 
and  is  printed  out  in  tabular  form  by  AVERAGE,  subroutine  PLOT  plots  out  the  data 
from  the  one  run. 
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VII.  RESULTS 


Figure  12  is  a  sample  of  the  output  of  the  chart  recorder.  This  example  is 
from  one  of  the  Iceland  boresight  runs.  The  patterns  were  recorded  at  11.  401 ,  30.  0, 
13. 996,  and  19.  0  MHz.  The  calibration  curves  for  each  frequency  are  shown  at  the 
end  of  the  run.  Since  the  data  are  smooth  (no  peaks  and  nulls),  the  readings  were 
taken  every  10  nmi.  In  Figure  13  there  is  a  sample  of  Starr  Hill  data  at  frequencies 
of  11.18,  15. 047,  17.  983,  and  25. 335.  These  curves  consist  of  several  peaks  and 
nulls,  and  for  this  example,  these  points  were  read  instead  of  the  10  nmi  points.  The 
data  read  from  these  plots  are  placed  on  data  sheets  (Figure  14).  For  each  frequency 
on  each  run  there  is  a  column  for  ground  distance  (NM),  divisions  on  the  chart  (DIV), 
and  attenuation  in  the  circuit  (dB).  There  is  also  at  least  one  calibration  curve  for 
each  frequency,  hi  most  cases  calibration  curves  were  run  both  before  and  after  the 
run. 


These  data  are  punched  on  computer  cards  and  serve  as  input  to  the  ANTPAT 
program.  This  program,  which  was  discussed  in  Section  VI,  has  two  outputs:  (1)  a 
printout,  and  (2)  a  data  tape.  A  sample  of  the  printout  (Figure  15)  lists  the  elevation 
angle  (degrees),  the  range  (nmi),  the  power  at  the  receiver  (dB),  attenuation  in  the 
circuit  (dB),  the  received  signal  (dB),  the  free  space  loss  (dB),  antenna  gain  (dB), 
and  the  loop  gain  (dB)  in  columns  1  through  8  respectively.  In  this  and  all  of  the  fig¬ 
ures  of  this  section  the  loop  gain  has  been  set  to  zero  (see  Section  V)  and  therefore 
the  data  represents  only  relative  gain.  The  magnetic  tape  output  was  used  for  further 
analysis  of  the  data. 

A  typical  plot  from  Program  TAPPLT  of  data  from  the  Panama  antenna  is  shown 
in  Figure  16.  The  two  curves  are  plots  of  different  runs  on  the  same  frequency  and 
the  same  azimuth  and  hence  represent  data  from  two  printouts  like  the  above.  The 
number  after  the  identifying  mark  is  the  date  of  the  run.  The  flight  plan  called  for 
data  runs  on  boresight,  ±5  degrees  from  boresight  and  ±10  degrees  from  boresight 
with  repeats  on  the  boresight  and  ±10  degree  runs.  This  plan  was  followed  almost 
completely  in  Panama  flights. 

The  average  and  standard  deviation  of  all  runs  on  the  same  azimuth  and  fre¬ 
quency  were  calculated  using  the  DATAVE  Program.  The  outputs  of  this  program 
include  both  a  printout  and  a  plot.  The  printout  (Figure  17)  lists  the  average  (dB), 
elevation  angle  (deg),  standard  deviation  (dB),  and  the  number  of  points.  The  number 
of  point"  is  given  so  that  we  may  assign  some  confidence  to  the  value  of  the  standard 
deviation.  Figure  18  shows  the  average  of  the  curves  in  Figure  16.  Average  curves 
of  the  Panama  antenna  were  plotted  for  each  frequency  and  heading.  Figures  19 
through  26  show  the  averages  on  boresight  and  ±10  degrees  from  boresight.  The  ±5 
degree  azimuths  are  omitted  for  clarity  of  presentation.  It  will  be  noticed  that  the 
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Figure  12.  Chart  Recorder  Output  - 
Iceland 


^  31/32 


iiiiiiiiiniiiiliiiiiiiiiiiii 


i  i  ii  fiiH  H  Hi  Hi  1 1 1 1 1 1  i  i  i  I H  i  tl  ii  ■  1 1 1 1  g  i  i  ■  J  g  i  1 1 U  lI  i  li  mi.1 


Jliltlft 


iiiniii'iiiiiiHiaiHiiiiiiiiniiimiHlWll  —1—iMBlMni 


SHUSH  r\STRlJIM  HVTH 


iiiiUri 


f 

o 


i 


I 


I  ,  , 

1  :  ! 

i 


i 


i 

i 

i 


ZOOOOOOOOOOOl 


0000000000900 


CM 


Ck. 

o 

o 


X  *<vcv.ocv^»*  o  If  o  o  c  c  o 

t-0«400»NrjlON<)'OniADlBNrl(0/\ntino»«o 

tfl  — 

K  2«0>lftnNri»0tl0'HOOfflJ>»f»r0r0W  I  •  CM  CM  fO  "> 


at  JtnNfllMNIMHHrlHHrll  •  »  »  •  I  I  I  t  I  ■  I  H 

Ml  K  I  I  I  I  I  I  I  I  I  I  I  I  • 


O  i/i 

a 


CM 


<3 

A 


•o 

7)  IN 

20 

O- 

rv  to 

-o  -o 

•O  ?V  GO 

to  ♦ 

x  in 

*o 

N. 

♦ 

rs. 

NO 

r-t 

to 

r4 

CM  r-t  > 

X 

UJ 

u»»Nhn»«« 

M  33  CM  -©  O 

cm  « 

Mt 

04 

O  fO 

O 

40  to 

K)  >£)  CM  CM 

X 

o  * 

• 

» 

• 

•  • 

t»  • 

•  •  • 

•  • 

• 

• 

• 

• 

• 

• 

• 

•» 

• 

• 

•  • 

• 

• 

• 

_» ♦ 

* 

tfOnNNHrlOO»lO 

co  x  in 

CM  rH 

N^OCM 

aH  O 

X 

** 

o 

t- 

r4 

H 

*4 

r< 

rt  H 

H  rt 

rtrtrt 

o  o 

o 

o 

o 

o 

O 

o 

O 

o 

a*- 

O  X 

X 

X 

z 

LU  rH 

O 

H 

H 

H 

rt  rl 

H  H 

H  H  H 

rH  rH 

H 

rt 

H 

rt 

H 

rt 

z 

< 

3 

Q. 

QC 

o 

x 

H 

c 

o 

o 

JlOO 

O 

> 

A  « 

A")  > 

H  O 

o 

rv 

V 

rl 

o 

o 

rvi 

r-t 

X 

> 

H  H  > 

V. 

«  rj 

It 

O 

X  wi 

C  Tt 

C3  O  W 

<  ir 

IT 

M' 

« 

IT 

< 

to 

r- 

« 

r- 

LU 

z  * 

• 

• 

• 

m  • 

•  • 

•  •  • 

•  • 

■ 

• 

•ft 

• 

• 

• 

■ 

• 

• 

• 

•  • 

• 

• 

»- 

C3  M 

CV 

r*> 

o 

CO  PN  l#N  CM 

»N«IOW 

O  W  fO 

rt 

O 

X 

nO 

to  rv  >o 

>C  Xl  IfN 

o 

•< 

-OlrfHHOO 

O  O 

O  O  »  O  O 

O 

OO 

CO 

CO 

®  N  rs  N 

o  -o 

0-0-0 

fv 

C9 

a 

/)  ri 

H 

rt 

rH 

H  H 

rH  rH 

z 

r-« 

a 

LU 

2  O 

CD 

o 

O 

o  o 

O  3 

0  0  3 

o  o 

o 

o 

o 

CD 

O 

o 

o 

o 

o 

O 

O  3 

o 

3 

X 

►- 

LU  ■ 

1 

• 

V 

1  1 

i  • 

o  c  o 

o  o 

CD 

c 

o 

O 

c 

o 

o 

o 

c 

o 

c  o 

a 

C 

LU 

>- 

•  •  • 

•  • 

• 

• 

• 

r 

» 

• 

• 

• 

• 

m 

*  • 

• 

• 

-U 

U 

•— 

t\  s\  n 

3  o 

XN 

n 

o 

3 

rv 

o 

3 

rv 

rv 

3 

ITS 

u. 

< 

r-t  rt 

rH  CM  CM  CM 

cm  cm  ro  ro 

'O  to 

ro  * 

X 

to 

tr 

X 

t— 

o 

o 

IO 

x>  fo 

o 

■c 

:> 

7  rs  ifi  * 

</>  *o  > 

rH  X 

o 

rv 

H 

X) 

o 

rvi 

r-t 

X 

> 

H  r-t 

t> 

Pv 

u 

C  K 

X 

w-< 

c 

C  r- 

C  c  r- 

x 

IT 

a 

X 

p^- 

K 

h  X 

rH 

* 

Q 

• 

• 

• 

• 

•  • 

•  » 

•  •  • 

•  • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

•  • 

• 

• 

o 

K> 

c? 

oo  rv 

IT  CM 

tr  a  rt 

r)  mh 

X 

tO 

rH 

rr 

X 

ro 

H 

M-t  i> 

rv 

IT 

X 

I— 

CM 

H 

r4 

H 

o  o 

o  o 

o  o  o 

O  3 

CD 

o 

CD 

O 

o 

o 

o 

o 

o 

O 

CD  3 

o 

3 

c 

► 

«-t 

rr* 

«-» 

v'  V* 

w-  r- 

■»t  v  rrt 

rt 

V- 

v* 

w* 

r- 

■ 

W 

v" 

r1 

r< 

^  mt* 

rH 

V 

o 

o 

o 

• 

< 

o 

LU  O 

C2* 

3 

o  o 

o  o 

o  o  o 

O  O 

CD 

o 

O 

O 

O 

O 

o 

o 

O 

3 

3  3 

O 

O 

ro 

CD  O 

Z  • 

O 

• 

C 

• 

o 

• 

o  o 
•  • 

o  o 
•  • 

o  o  o 
•  •  • 

o  o 

•  • 

CD 

• 

o 

■ 

O 

• 

c 

• 

3 

• 

O 

• 

o 

• 

o 

• 

o 

• 

O 

• 

O  O 

•  • 

o 

• 

O 

• 

•<OCOi-IOOOOOOOOCO»-tOOOOOOC'Opsljr'CVlOir. 


<r 

a: 

PO  CM  r-t  O  O 

x  rv 

-c  in  r 

ro  cm 

r-l 

X 

rv 

«o  rv 

to 

CM  r-t 

rH  rH 

rH 

tl 

CM 

t\J  'M  CM  CM  r-l 

H  rt 

r-t  r-t 

r-t 

r-t  -4 

r4 

u 

X 

z 

LU 

IU 

> 

a 

m  rv 

-O  •O  lfN  O  OO 

rrt  X 

CM  O 

X 

(MON®  O 

o 

CM  CM  50 

♦ 

r-t  ro 

rv 

'O  M- 

o 

Uj 

-j 

rt 

«=  C  r-l  CM  KD  IT  <> 

tc  c 

r-t 

r:  x 

X  IT 

c 

ITv 

ro  to 

rv 

r-t  rv  c 

C  CM 

rH  X 

LU 

o 

a 

• 

•  • 

•  • 

• 

•  • 

•  • 

• 

■ 

r  • 

• 

• 

•  m 

•  m 

•  • 

oc 

LU 

z 

« 

I 

r-t 

r-t 

r4  rH 

r-t  CM 

ro 

to 

▼  LTV 

o 

to  X 

X  CM 

x  "<r 

u. 

o: 

< 

r-t  r-t 

rH  CM  CM  M- 

36 


I 


Figure  15.  Program  ANT  PAT  -  Tabular  Output  -  Panama 


-24  .O 


rAN-V 


28.0  -24.0 


-20.0  -16.0 

□B 


-I  2 .0  -8.0  -4.0  .O 


Figure  16  -  Program  TAPPLT  - 
Plot  Output  (Panama) 


37/38 


I 


FREQUENCY*  30 
AVE  GAIN 
-36,0001 
-29,6162 
-25,6206 
-23,1000 
-22,115$ 
-20,594$ 
-19,3517 
-17.4197 
-15,4039 
-14,3010 
-13,1356 
-11,3356 
•9,6034 
•0,6331 
•5,9066 
•4,7304 

•3,0700 

•3,7048 

•3,1479 

•2,5084 

•3,2464 

-.5137 

•1.0701 

•2.3234 

•2,8997 


POLARIZATION  no, 
ANGLE 
-.1653 
-.0646 
.0291 
,1400 
,2600 
.3709 
.5111 
,6620 
.0220 
,9955 
1,1026 
1,3920 
1.6322 
1,0737 
2,5034 
3,0292 
3,5929 
4,3243 
5,3211 
6,7769 
9,1423 
13,7134 
16,0317 
10,0703 
22,2359 


5  ANGLE*  6 

STD, DEV, 
n 
0 
0 

.104? 

,1003 

.3327 

.0610 

.5776 

1.2309 

.5810 

.550? 

,7l49 

.2700 

.1024 

.002? 

,077§ 

.3423 

.414? 

.631? 

.0709 

,0536 

.0603 

.5221 

,2301 


POINTS 


Figure  17.  Program  -  DATAVE  -  Output  (Panama) 
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-10  degrees  from  boresight  (356°)  Is  consistently  higher  than  the  boresight  (006°). 
There  is  a  trend  of  increased  pattern  gain  to  the  left  of  boresight.  (This  trend  is  borne 
out  by  the  ±5°  azimuths  also. )  Thus  the  direction  of  the  main  beam  of  the  pattern  dif¬ 
fers  by  some  10  degrees  from  the  direction  in  which  the  antenna  is  pointing.  The 
main  reason  for  this  is  the  proximity  of  sea  water  to  the  left  side  of  the  antenna. 

The  next  step  was  the  reduction  of  these  patterns  to  a  form  that  can  be  used  in 
the  mode  loss  computer  program.  There  are  a  number  of  methods  that  could  be  em¬ 
ployed  to  do  this.  At  first,  consideration  was  given  to  fitting  these  curves  with  poly¬ 
nomial  functions  of  the  elevation  angles.  If  the  data  could  be  matched  with  a  set  of 
low  order  polynomials  it  would  mean  a  savings  in  computer  time  and  space.  This 
method  turned  out,  however,  to  be  less  efficientthan  a  simple  table  look-up.  This  latter 
method  consists  of  (1)  setting  up  a  set  of  tables  of  gain  versus  elevation  angle ,  (2)  using  this 
table  to  interpolate  to  the  elevation  angle  of  interest  on  the  two  frequencies  adjacent  to  the 
frequency  of  interest,  and  (3)  interpolating  between  the  two  values  thus  obtained.  Therefore 
the  gain  can  be  calculated  at  any  elevation  angle  and  frequency  in  the  range  of  the  data . 

Sample  results  of  the  measurements  of  the  Thule  antenna  are  shown  in  Figures 
27,  28,  and  29.  These  results  are  for  typical  low,  mid,  and  high  frequencies.  These 
results  are  completely  at  variance  with  the  results  obtained  on  the  Panama  antenna. 

The  Thule  antenna  siting  is  the  cause  of  this  dissimilarity.  This  antenna  is  located 
about  one-half  mile  above  sea  level  (see  Figure  3).  The  lobes  in  this  pattern  are 
caused  by  (1)  interference  between  the  antenna  and  its  image,  and  (2)  diffraction 
around  the  sharp  ridges  in  front  of  the  antenna.  A  discussion  of  this  problem  may 
be  found  in  Reference  1.  Figures  30  through  37  are  the  average  patterns  of  the  Thule 
antenna  on  boresight.  The  patterns  at  the  other  azimuths  are  similar.  Looking  at  the 
patterns  at  different  frequencies,  it  can  be  seen  that  the  peaks  and  nulls  occur  at  dif¬ 
ferent  angles  in  each  case.  Because  of  this ,  we  cannot  directly  interpolate  between 
frequencies.  Instead,  a  method  was  derived  which  fits  the  data  in  terms  of  the  num¬ 
ber  of  the  null. 10  Using  this  method  an  interpolation  scheme  was  established  which 
takes  into  account  the  location  of  these  peaks  and  nulls. 

Typical  data  curves  for  the  Iceland  antenna  patterns  are  shown  in  Figure  38  and 
their  average  is  in  Figure  39.  Although  the  flight  plan  called  for  the  same  number 
runs  as  at  other  sites,  less  valid  data  runs  were  obtained.  Several  data  runs  could 
not  be  used  because  of  instabilities  of  the  voltage  levels  in  the  recorder.  The  average 
curves  for  the  Iceland  antenna  are  shown  in  Figures  40  through  47.  The  curves  are 
shown  for  boresight  (260°)  and  ±10°  off  boresight.  There  is  a  tendency  for  the  +10° 
bearing  (270°)  to  have  a  higher  gain  than  the  boresight  and  for  both  of  these  to  be 
higher  than  the  -10°  (250°)  bearing.  Since  the  gains  for  the  260°  and  270°  azimuths 
are  close  together  it  appears  that  the  peak  of  the  beam  is  offset  less  than  10  degrees 
to  the  north  from  the  boresight  of  the  antenna.  As  in  the  case  of  the  Panama  antenna, 
the  ocean  is  closer  on  the  side  to  which  the  beam  is  offset. 
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Figure  29.  Program  TAPPLT  -  Hot 
Output  (Thule  -  30  MHz) 
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Figure  30.  Thule  Antenna 
Patterns  -  9  MHz 
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Figure  36.  Thule  Antenna 
Patterns  -  26  MHz 
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Figure  40.  Iceland  Antenna 
Patterns  -  9  MHz 
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Figure  41.  Iceland  Antenna 
Patterns  -  11  MHz 


89/90 


DB 


Figure  42.  Iceland  Antenna 
Patterns  -  13  MHz 
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Patterns  -  17  MHz 
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Figure  44.  Iceland  Antenna 
Patterns  -  19  MHz 
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It  was  discovered  that  the  Iceland  antenna  pattern  data  was  subject  to  tidal  ef¬ 
fects.  Figure  48  shows  two  data  runs  flown  on  the  same  day,  the  first  at  high  tide  and 
the  second  at  low  tide.  The  difference  in  water  level  between  these  two  cases  is  ap¬ 
proximately  eight  feet.  These  variances  in  the  pattern  are  too  large  to  be  disre¬ 
garded.  Before  utilizing  these  data,  a  method  must  be  derived  accounting  for  these 
effects.  The  simplest  method  would  be  to  obtain  an  average  gain  and  then  consider 
the  tide  effects  as  a  perturbation  on  this  average. 
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VIII.  FUTURE  CONSIDERATIONS 


The  gain  measurements  of  the  antennas  at  the  Starr  Hill  receive  site  have  not 
been  completed.  The  measurement  technique  discussed  herein  will  be  employed,  and 
it  is  intended  that  this  phase  of  the  program  be  completed  during  the  spring  and  sum¬ 
mer  of  1968.  In  addition,  the  calibration  of  the  aircraft  loop  antenna,  required  to 
obtain  absolute  gain,  will  be  accomplished  during  this  same  period. 

A  calibration  oi  the  azimuth  of  arrival  measurement  system  (antennas  and 
equipment)  located  at  Starr  Hill  is  tentatively  scheduled  during  the  summer  and  fall 
of  1968.  Transmissions  from  a  Convair  C-131  aircraft  will  be  used  in  these  tests. 
The  results  of  the  Starr  Hill  measurements  will  be  reported  on  in  the  future. 
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EXPANDED  LITTLE  IDA  PROGRAM 

The  Expanded  Little  IDA  effort  involves  the  collection  and  analysis  of  high  fre¬ 
quency  propagation  data  that  will  assist  in  the  design  and  operation  of  over-the- 
horizon  detection  systems.  Several  sophisticated  propagation  experiments  are  being 
conducted  over  three  2000-nautical  mile  paths.  The  paths  are  (1)  a  transauroral  path 
between  "P"  Mountain,  Thule,  Greenland,  and  Starr  Hill  Test  Annex,  (2)  a  "parallel- 
to-the-auroral  zone"  path  from  Site  Dye  5,  Iceland  to  Starr  Hill,  and  (3)  a  temperate/ 
semitropical  path  between  Coco  Solo  Test  Annex,  Panama  Canal  Zone  and  Starr  Hill. 
The  Thule,  Iceland,  and  Panama  sites  transmit  their  signals  to  the  Starr  Hill  central 
receiving  site. 

The  experiments  and  investigations  that  are  being  performed  include: 

A.  1.  Mode  Reliability 

The  object  of  this  experiment  is  to  determine,  (1)  the  availability  of  the  lower 
order  dominant  modes  (one  hop  F2,  2  hop  F,  etc. ),  (2)  the  maximum  and  lowest  ob¬ 
served  frequencies  used  by  the  individual  propagating  modes;  and  (3)  the  "available 
frequency  aperture"  or  range  of  usable  frequencies  for  each  mode  of  propagation. 
These  parameters  are  being  measured  as  a  function  of  time  of  day,  month,  season, 
and  path.  This  experiment  is  designed  to  answer  the  following  questions:  "What 
modes  are  available  and  are  reliable  over  the  regions  of  interest  at  any  given  time?", 
"At  what  frequency  or  band  of  frequencies  should  a  radar  operate  to  illuminate  a 
specified  region?",  and  'What  upper  and  lower  frequencies,  or  what  range  of  oper¬ 
ating  frequencies,  should  be  specified  for  the  design  of  a  specific  system?". 

A.  2.  Mode  Loss 

The  goal  of  the  mode  loss  experiment  is  to  determine  the  losses  associated  with 
the  individual  modes  that  propagate  over  the  three  paths.  Certain  modes,  with  high 
availability,  may  be  undesirable  due  to  the  high  losses  associated  with  them.  Be¬ 
cause  of  signal -to-noise  considerations,  there  is  a  need  to  ascertain  the  loss  suffered 
by  each  ;ode  and  the  availability  of  the  mode  as  a  function  of  time  of  day,  month, 
and  season.  In  addition  to  the  loss  measurement,  an  effort  is  being  made  to  separate 
and  identify  the  predominant  sources  of  loss,  such  as  polarization  effects,  polar  cap 
and  auroral  zone  absorption,  focusing,  and  Es  layer  masking. 

The  analysis  of  propagation  loss  due  to  various  mechanisms  requires  the  de¬ 
termination  of  the  total  loss  between  transmitter  and  receiver  and  the  separation  of 
known  losses  and  gains. 
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This  is  further  demonstrated  by  considering  the  following  equation: 


L. .  +  L  -G 
of  P  P 


where 

Hot  ”  *s  total  system  loss  between  the  transmitter  and  receiver.  This 

quantity  is  determined  from  the  measurements  and  is  computed  directly 
from  the  RF  power  delivered  to  the  transmitting  antenna  line  and  the 
RF  power  delivered  to  the  receiver. 

-  is  the  total  antenna  power  gain  for  a  given  mode  and  is  equal  to  Gp  =  10 
log  gr  +  10  log  gt  where  gr  is  the  receiving  antenna  power  gain  for  a 
particular  mode,  and  gt  is  the  transmitting  antenna  power  gain  for  a 
particular  mode.  The  quantity  Gp  is  determined  from  the  calibrated 
antenna  tests. 

-  is  the  basic  free  space  loss  or  the  loss  due  to  the  spreading  of  the 
electromagnetic  wave  in  space  and  is  defined  as: 

_  4»a2 

h>f=  x 

where  R  -  equivalent  distance  between  transmitter  and  receiver,  and 
X  =  RF  wavelength. 

-  is  termed  the  propagation  loss  and  includes  the  ionospheric  losses  and 
ground  reflection  losses.  The  value  of  Lp  and  its  variability  on  an  in¬ 
dividual  mode  basis  as  a  function  of  path,  time  of  day,  month,  and 
season  is  the  ultimate  objective  of  this  experiment. 

It  is  therefore  evident  that  the  antennas  and  transmit/receive  systems  must  be 
well  calibrated  so  that  the  system  losses  and  antenna  gain  parameters  may  be  re¬ 
moved  from  the  measured  total  loss.  In  particular,  the  gain  patterns  must  be  well 
known  as  a  function  of  elevation  and  azimuth  angles. 

This  experiment  is  being  conducted  over  paths  that  have  a  2000  nm  range  wmch 
is  the  limiting  distance  for  propagation  via  the  one-hop  F  (1F2)  mode.  When  the  1F2 
mode  is  propagated  over  these  ranges,  it  is  received  at  extremely  low  elevation 
angles  (from  -1*  to  +4°).  On  the  other  hand,  the  two  hop  F2  mode  is  received  at  ele¬ 
vation  angles  between  12*  and  20*.  Therefore,  the  antenna  gain  patterns  must  be 
known  as  a  function  of  elevation  angles  that  vary  from  -1°  to  +20°.  In  addition,  the 
elevation  gain  pattern  must  be  known  as  a  function  of  azimuthal  angle  since  the  indi¬ 
vidual  modes  may  be  propagated  via  off-great-clrcle  paths. 
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A.  3. 


Polarization 


The  objectives  of  this  experiment  are  (1)  to  acquire  fading  statistics  (cross¬ 
correlation)  between  signals  received  on  cross-polarized  antennas  so  as  to  determine 
whether  or  not  the  losses  associated  with  the  individual  modes  are  polarization  sensi¬ 
tive,  and  (2)  to  determine  whether  or  not  any  of  the  propagation  paths  exhibit  a  pre¬ 
ferred  polarization  as  a  result  of  magneto-ionic  effects  of  the  ionosphere  on  wave 
propagation.  To  conduct  the  measurement  program  for  this  experiment,  the  receive 
site  has  been  instrumented  with  two  arrays,  (oriented  North  and  South)  consisting  of 
independent  log-periodic  vertically  and  horizontally  polarized  antennas.  To  ade¬ 
quately  perform  the  measurements,  each  antenna  oeam  should  exhibit  a  reasonable 
cross-polarized  response  of  -20  db  or  less.  In  addition,  it  is  required  that  the  two 
independent  orthogonally -polarized  antenna  beams  be  matched,  to  within  a  few  db, 
over  the  expected  vertical  arrival  angles  (0°  to  20°). 

A.  4.  Azimuth  of  Arrival 

The  purpose  of  this  investigation  is  to  determine  the  azimuth  of  arrival  char¬ 
acteristics  of  the  dominant  lower  modes  of  propagation.  The  extent  and  variability 
of  off-great-circle  path  propagation  of  the  individual  modes  is  to  be  determined.  At 
the  Starr  Hill  receive  site  each  antenna  array  contains  two  log-periodic  vertically 
polarized  dipole  antennas  (4-64  Me) ,  arranged  at  a  constant  spacing  of  two-thirds  of 
a  wavelength  for  use  as  an  interferometer.  The  azimuth  of  arrival  measurement  sys¬ 
tem  (antennas  and  receiver  system)  must  be  well  calibrated.  The  results  of  this 
calibration  (which  will  also  employ  aircraft  measurement  techniques)  will  be  reported 
on  at  a  later  date. 

A.  5.  Noise  and  Interference 

This  effort  requires  determining  as  precisely  as  possible  the  noise  and  inter¬ 
ference  environment  within  which  an  OHD  system  would  operate.  The  investigation 
involves  two  separate  measurement  techniques. 

The  first  technique  involves  the  collection  of  omni-directional  noise  data  at 
the  Starr  Hill  site  by  the  standard  CCIR  (Consultative  Com.  littee  on  International 
Radio)  measurement  approach  which  is  used  in  the  world-wide  network  of  noise  Re¬ 
cording  stations.  This  method  utilizes  a  22  ft.  whip  antenna  in  conjunction  with  an 
ARN-2  ESSA  (Environmental  Science  Services  Administration)  receiving  and  record¬ 
ing  system  which  measures  noise  level  at  2.  5,  5,  10,  and  20  MHz.  In  addition, 
directional  noise  data  is  also  measured  using  a  separate  ARN-2  receiving  system 
and  the  log-periodic  vertical  and  horizontal  antennas.  The  measured  omni¬ 
directional  CCIR  noise  data  is  compared  with  (1)  predictions  published  by  ESSA  and 
(2)  the  directional  CCIR  type  noise  data  to  determine  if  the  noise  is  sensitive  to  pol¬ 
arization,  antenna  directivity,  and  also  to  obtain  a  better  understanding  of  the  effects 
of  antenna  characteristics  (high  angle  reception  vs.  low  angle  reception)  on  noise  and 
interference  data. 


The  second  noise  measurement  technique  uses  a  special  purpose  dual  channel, 
high  dynamic  range  <120  dB)  noise  receiver  with  a  frequency  range  of  4  to  40  MHz. 
This  receiver  has  tins  capability  of  scanning  250,  four  KHz  steps  in  any  of  the  36,  one 
ICHz  bands.  This  equipment  has  selectable  integration  times  of  0. 3,  3. 0  and  500 
secs.  The  500-second  Integration  time  is  also  common  to  the  ARN-2  receiver.  The 
parameters  of  this  receiver  are  comparable  to  those  presently  under  consideration 
for  use  in  proposed  systems.  Noise  data  as  a  function  of  frequency,  polarization, 
directivity,  and  bandwidth  are  being  collected  with  this  equipment  at  the  Starr  Hill 
Site.  This  type  data  is  more  meaningful  than  the  COIR  type  for  application  to  OhD 
problems  since  the  data  is  gathered  on  the  entire  HF  band,  as  well  as  on  directional 
antennas. 

A.  6.  Spectrum 

The  objective  of  this  effort  is  to  determine  the  effect  of  the  ionosphere  and 
other  variables  on  signal  coherence  and  signal  doppler.  Signal  coherence  information 
pertaining  to  signals  propagating  over  each  of  the  paths  is  required  so  that  the  system 
designer  may  determine  the  optimum  processing  format.  That  is,  he  may  need  to 
utilize  coherent  Integration  to  achieve  a  required  signal-to-noise  ratio.  The  signal 
doppler  information  is  required  to  determine  the  system's  capability  in  detecting 
moving  targets  while  operating  in  a  doppler  processing  mode.  These  measurements 
are  made  on  an  individual  mode  basis  through  use  of  a  spectrum/analyzer  magnetic 
drum  processor. 

The  experiments  discussed  above  are  being  conducted  simultaneously  over  the 
three  paths  of  interest.  The  interdependence  of  the  various  parameters  is  also  under 
investigation  in  order  to  better  define  the  overall  effect  on  radar  signals  propagating 
over  the  regions  of  interest.  In  the  past,  these  experiments  have  been  performed  as 
separate  measurements  in  themselves.  Under  this  program,  defining  the  relation- 
shlpo  existing  between  the  parameters  measured  is  a  major  task.  For  example, 
noise  and  loss  have  been  treated  separately  and  independently  in  previous  endeavors. 
On  the  other  hand,  there  is  good  reason  to  believe  that  these  quantities  are  dependent 
and  are  negatively  correlated,  i.  e. ,  when  the  absorption  is  high,  the  noise  may  be 
low.  The  correlation  between  the  loss  and  noise  data  is  also  being  determined  u  a 
function  of  the  available  frequency  aperture  (range  of  usable  frequencies)  so  as  to  be 
more  meaningful  for  systems  application. 

From  the  discussion  above,  it  is  evident  that  the  systems  employed  in  the 
measurement  program  must  be  well  calibrated.  Since  the  antennas  are  a  major  sys¬ 
tem  component,  and  since  the  measurements  are  being  made  at  several  frequencies, 
a  detailed  antenna  gain  measurements  program  is  required  for  the  experimental  re¬ 
sults  to  be  meaningful. 

Detailed  descriptions  of  the  experiments  and  instrumentation  used  are  re)>orted 
on  in  References  1  through  6. 
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PROGRAM  ANTPATT 


PROGRAM  ANTPAT 

DIMENSION  C<50),  CC<50),  GR(IOO).  0(100),  AT(HO),  DB(100) 
DIMENSION  AA(9) *  AB(9) 

DIMENSION  NF<13) 

DATA! <NF<J), I *1,13 >*9.11, 13. 15. 16. 17. 19. 21. 22, 25, 26, 29, 30) 
DIMENSION  XX<100).YYtl00) 

TYPE  REAL  LGG10F 

CALIBRATION  DATA  (DBM) 

DATA! <CC< I), I *1,11) *115., 110. ,105., 100. ,9a, ,90. >85. ,80 . , 73. , 

1  70. .65.) 

DATA((AA(I),  1*1,9)*  -1.,  0..  1.,  2.,  3.,  4.,  6.,  10..  14.) 

DATA ( ( A0 ( I ) ,  1*1.9)*  .7,  .45.  .31,  .22.  .16,  .13,  .09,  ,05.  ,03) 
LOG10F (X) * . 43429*L0QF ( X ) 

REMIND  2 

STATEMENTS  407  THROUGH  12  ARE  TO  BE  REMOVED  FOR  THE  FIRST  RUN 
(WHEN  NO  DATA  IS  ON  TAPE  2) 

THIS  SECTION  LOCATES  END  OF  DATA  ON  TAPE  2 

407  READ  INPUT  TAPE  2 , 766, NRUM , NPOL . NDATE , NTHD, NRvN . NT  I  ME , NREC • ALTKF , 
1  FREQ. POME 
IF (NRUM  .EQ.  99)  12,75 

75  READ  INPUT  TAPE  2.777.NPT, (XXC I ) , YY( I ) » 1*1, NPT) 

GO  TO  407 
12  CONTINUE 
BACKSPACE  2 
NFREQ1*0 
NTHD1*0 
CNMKM*1. 83323 
CKFNMaO. 164462 


SET  POLARIZATION - 

TRANSMIT 

RECEIVE 

--SITE 

1  * 

H 

w  ~ 

H 

--STARR  HILL 

2  * 

H 

m  m 

V 

--STARR  HILL 

3  * 

V 

H 

--STARR  HILL 

4  * 

V 

-- 

V 

--STARR  HILL 

5  • 

V 

-  - 

V 

--panama 

6  • 

V 

-  - 

H 

--PANAMA 

7  « 

V 

-- 

V 

—  THULE 

8  • 

V 

-- 

H 

—  THULE 

9  * 

V 

-  - 

V 

--ICELAND 

0  « 

V 

-- 

H 

--ICELAND 

NPOl*l 

RENM*3443 , 2 

RE2*RENM**2 

FPI*4. *3.14159 

CRADDEG*57.2938 

CONSENT  LOSS  FACTOR 
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CLOSSa60.-2Q.*LOGl£)F<  30  0.  >  *20  .  *LOG10F  ( FP I  > 

READ  FIRST  DATA  CARD 

76  READ  77. NDATE, NTHD. NRUN.NTIME, ALTKF 

77  FORMAT! 16. 13.12. 14. F10.0) 

IF(NDATE)  10,10,101 

10  NRUM*99 

WRITE  DUMMY  RECORD  ON  TAPE  2  ,  REWIND  THE  rAPE  ,  END  THE  PROGRAM 

766  FORMAT! 12. 1 1. 16, 13, 12. 14, H.3E16.9) 

WRITE  OUTPUT  TAPE  2 , 766 , NRUM , NPOL , NDATE , NTHD . NRUN , NT  I  ME , NRfcC . 

1  ALTKF.FREO.POWE 
6  REWIND  2 
STOP 

777  FORMAT!  I2,6E16.9/(6E16.9)  ) 

101  nai_t«altkf 
IF<NALT)  102.102,103 

102  ALTKF *30.0 

READ  SECOND  DATA  CARD- <  NUMBER  OF  DATA  PO I  NTS , RECE I VER  NO ., FREQUENCY  ) 

103  READ  ll.NPOINT.NREC.FREQ 

11  F ORMAT ( 2  I  2 , F10 . 0  ) 

NFREQ*FREQ 

IF(NFREQ.EQ.NFRE01. AND. NTHD. E0.NTHD1 )  701.702 

701  NRUM»0 

GO  TO  106 

702  NRUU«1 
NFREOlxNFREQ 
NTHDlsNTHP 

108  CONTINUE 

READ  THIRD  DATA  CARD  - (POWER  IN  DBM) 

READ  88.  POWE 
88  FORMAT ( FiO • 2 ) 

READ  CALIBRATION  DATA 

104  READ  5>5 ,  ( C (  I  )  ,  I«*,ll) 

55  F ORMAT ( 16F4 , 1 ) 

PRINT  DATA  HEADING 

PRINT  22.FREQ.NTHD.NRUN.NREC. AlTKF .NDATE, NT  I  ME, POWE, NPOL 
22  FORMAT ( 1H1 , // . 5X , lOHFREQUENCYs , F 10 . 4 , 3HMCS , 5X » 12HTRUE  HEADING, 15, 
C5X.7HRUN  NO. , I  3 , / , 5X . 8HRECE I VER , J 3 , 5X . SHALT  I TUDE  , F5 . 1 , 5HKFEET , 5x , 
C4HDATE. 18.4X.4HTIME, 1 8 , 4X , 6HP0WER* . FIO . 2 , 2HDB , 7x , 4HP0L= « 14) 

PRINT  66 
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6*  FORMAT ( // 1 5X  » 5HANGLE  *  5X . 5HRANGE .8X.2HDB.5X.5HATTEN.4X. 6HS I GNAL , 
C5X.10HSPACE  LOSS .5X.15HS1 GNAL  STRENGTH. IX. 9HL00P  GAIN) 

READ  —GROUND  RANGES— DIVISIONS— ATTENUATIONS 

READ  33»<GR(I).  I«1.NP0INT) 

READ  33. (D( I ) i  Isl.NPOINT) 

READ  33. ( AT(  I ) .  I»1,NP0INT) 

33  FORMAT (15F5.2) 

NPOINTlsO 

I«1 

NNFREOaFREO 

SET  NFF  TO  EOUAL  THE  NUMBER  OF  THE  FREQUENCY* (FROM  1  TO  8) 

DO  301  Nisi. 13 

IF (NNFREQ  .E0.  NF(ND)  302.301 

301  CONTINUE 

302  NFFsNI 

THIS  SECTION  (7HR0UGH  STATEMENT  906.)  IS  REPEATED  FOR  EACH  DATA  POINT 
DO  906  Jsl.NPOINT 
CONVERT  DIVISIONS  TO  DB'S 

IF(D(J)-C(1)>  7.2.6 

6  DO  7  Isl.ll 

IF (D( J) -C ( I ) )  3.2.7 

7  CONTINUE 
PRINT  555 

555  FORMAT ( 19H  OUT  OF  CALIBRATION) 

GC  TO  906 

2  DB(J)sCC(I) 

GO  TO  106 

3  XtD(J)-C( 1-1) 

V  s  C(I)  -  C(I-l) 

DB(J)sCC( I -1 > ♦X*CCCC I)-CC< I-l))/Y 

CALCULATE  RANGE  AND  TAKECFF  ANGLE 

106  GRANNMsQR(J) 

DBRECsDB(J) 

ATTENsAT(J) 

ALTNM«ALTKF*CKFNM 

AH  s  RENM  *  ALTNM 

AH2  s  AH««2 

DRAD  «  GraNNM/RENM 

R2»RE2*AH2“2.#RENM*AH*C0SF(DRAD) 

SRANNM  S  SQRTF (R2) 

srankm  s  srannm*cnmkm 

A2  s  ALTNM*s2 
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SINTOA  •  <2.«RENM*ALTNM*A2-R2)/(2.*RENM«SRANNM) 
TOAR  «  ASINF (SINTOA) 

T0AxT0AR*CR ADOEG 


CORRECT  TAKEOFF 


ANGLE  FOR  RF^RACT ION 


401 

402 

403 

404 

405 


IF ( TO A  . GT  >  14  )  405, 

DO  404  10  s  1,9 

1 F ( TOA  -  A A ( 1 0 ) )  402 

JO  x  IQ  -  1 

TOA  x  T0A*AB(  I Q ) - ( AS 

GO  TO  405 

TOA  x  TOA  ♦  AB ( 10) 

GO  TO  405 

CONTINUE 

CONTINUE 


401 

403,  404 

I Q ) " AB ( JQ ) )*(AA( ]Q)-T0A)/(AA( I Q ) - A A ( JO  )  ) 


CALCULATE  signal  STRENGTH 


SUBROUTINE  LOOP  CALCULATES  GAIN  OF  LOOP  RECEIVING  ANTENNA 
(THIS  ROUTINE  IS  NOT  USED  AT  PRESENT) 

CALL  LOOP (NFF ,T0A,22) 

FSLxCLOSS*20.*L0G10F(FREQ)*20.*LOG10F(SRANKM) 

TSIGxDBREC-ATTEN 

SSTxFSL-TSIG-POWE*ZZ 

PRINT  DATA 

PRINT  44, TOA.GRANNM.DBREC, ATTEN , TS 1 G  ,  FSL , SST , ZZ 
44  FORM4T(5F10.2.5X,F10.2,10X,2F10.2> 

RECORD  DATA  POINT  IN  XX  AND  Y  Y  ARRAYS 

I F ( TO  ^  . GT . 25 . )  906,  705 
705  NP0INfl»NPOlNTl+l 
XX(NPOINTl)xSST 
YY(NPOINT1)3TOA 
906  CONTINUE 


WRITE  DATA  OUTPUTS  ON  TAPE  2 

WRITE  OUTPUT  TAPE  2 . 766 . NRUM , NPOL . ND ATE , NTHD , NRUN , NT  I  ME , NR6C > 
1  Al.TKF  ,  FREQ  ,  POWE 

WRITE  OUTPUT  TAPE  2 , 777  , NPO I NT1 , < XX (  I  )  , YY ( I ) » J *1 , NPO I NT1 ) 

READY  FOR  NEXT  DATA  RUN 

GO  TO  76 
END 

SUBROUTINE  LOOP(NFF,  TOA,  ZZ) 

ZZxO.O 
12  RETURN 
END 
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PROGRAM  TAPPLT 


PROGRAM  TAPPLT 

DIMENSION  mlOO)?  YYliADJ _ 

DIMENSION  NF(13),N*N(J(35),NPU(iO> 

I3ATAKNF (lli 1*1. 13)  *9 .11j13i15,’16  >  17>l9, 21. 22 » 25, 26. 29 . 30 ) 

C  STARR  HILL 

DATA<  < N A N fi <  1 )  i  I*l,7)i357i2, ? 1 3*. 183j  180,191) 

DATA ( ( NANO (I),  I«8.;>0)«352, 12*17,26, 31, 41»4*,l76,l98i 
1  917.910*919.920) 

C  MANAMA 

DATA { (NANG { I ) ,  I«2142l)  )i396>l»6.11.16) 

C  THULE 

DATAC (NANG( I ).I«26. 30)9179, 194.109. 194.199) 

C  ICELAND 

DATA(INANGU)?  1*31,3  J)42>0, 295, 260,265. 273  ) 

C 

DATA( (NPL( I ).l>1.10)«l42i3, 4 15,617.8,9,0) 

C 

.  DIMENSION  1FREQ(13),1HEAU(35),1PUL(10) 

DATA((IFREQ(I).I«1.'13)»2H9  ,  2HH ,  2M13, 2H15, 2H10 , 2H17 . 2H19 , 2H21 , 

C  2H22«2M?5,2H26,2H{9,2H3QJ 
C  STARR  HJLL 

DATA((IHSAD(I).I«1.7)«3H357,3H002.3H007,3H036,3H163,3H188.3H193) 

DATA((IHFAD(I).I*0V2O)«3H392,3HO12,3HO17,3HO26.3HO31,3HO41.3HO46. 
1  3H170. 3H190, 3H917 , 3H918, 3H919, 3H920 ) 
c  Panama 

DATA((IHEAD(I),I*2i,25)*3H356,3H001,3H006,3HOU.3H016) 

C  THULE 

DATA((IHEAD(I).I*26,3O)*3Hl7V.3H104,3H189,3H194,3H199) 

C  ICELAND 

DATA(  ( IHEO(  I ) ,  l«3l,  I5)«3H250 ,3H25b,3H260,3H265,3H270) 

C 

DATA((IPOL(l),I*l»lO)*5HSH-HHf5H6H-HV,5HSH-VH,5HSH-VV,5HPAN-V, 

1  5HPAN»X,5HTHU-V.5HTHU-X»5HlCEiV,5HICE-X) 

C 

YDAT«19,0 
SYHB  •  16, 

S^T  COUNTS/IN.,  BOARD  SCALESi  INTERPOLATION  DISTANCE 

CALL  PLTPT (300 • ,300 . .2) 

CALL  PLTPT (4, ,3. ,13) 

CALL  PLTPT(5, 0,0.0,14) 

START  MAIN  SEARCH 

1  READ  li.KK,K,K^K 
11  FORMAT  2312) 

IF  (K)  1110,U10,lll 
C  SET  PLOTTER  IN  57AN0SY  BEFORE  ENDING  PROGRAM 
1110  CALL  PLTPT(0, , 0 . ,5) 

STOP 

111  CONTINUE 
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REWIND  2 
NG0*JL  . 

911  READ  INPUT  TAPE  2 , 76A . NRUM , NPOL « NDATE , NTHD , NRUN , NT IMF. . NREC * ALTKF , 
l  FR6Q.P0WE 

766  FORMAT!  12, 11,  16, 13, 13, 1 4,12, 3El6. 9) 

777  FORMAT! 12, 6E16. 9/(*El<j ,  9  )  )  . 

NFR&QsFREQ 

I F !  NFREO ,  EQ ,  NF !  K_l-tiHD-.lslRQL.by .  N?V1JW )  .4NJ3.NTMD.tQ.  NANG !  KKK ) )  77,4 
4  I F ! NRUN , EQ  ,  99)  10,8 

77  READ  INPUT  TAPE  2 , 777 ,  NPT , < *X { I ) , m 1 1, ! •! . NPT > 

NPO I NT1*NPT 
1 F ( NGO . EO  .  1  )  706.7061 

SET  PLOTTER  IN  STANDRY 

706  CAUL  PLTPTIO.  -0.  ,5) 

NGQ*0 

IF  ( XX ! 1 )  -  XX!NP01NT1) )802,  801,  801 

801  RMIN  «  XX ( NPOINT1 ) 

GO  TO  803  . 

802  RMIN  «  XX  < 1 ) 

803  CONTINUE 

SET  PLOTTER  OFFSET 

NRM INaRMlN  —  7 . 0 
NXOFSTs-! (NRMIN-1) /4> 

XOFST«NXOFST 

CALL  PLTPT { XOFST , 1 , , 1 ) 

XMOFST*-XOFST 

FRMIN*XMOFST*4.0 

PLOT  SCALE 

CALL  PLTSC(FRMIN,4'f  ,10,,  1,-3,  ,3.  ,8.  .1) 

SELECT  LINE  MODE 

CALL  PLTPT (1 , , 0  . , 19 ) 

SELECT  SHORT  LINE 

CALL  PLTPTd.,0,,4) 

INTEG«4HHL-AD 

INTEG2«6HFREQ 

INTEG3«8HDEG 

INTEG4«6HDB 

INTEGX«8H- 

XSTARTsFRMIN-5,0 

PLOT  HEADING  INFORMATION 

CALL  PLTAL!XSTART,i2, ,0.25,INTGG3,3> 
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XSTART«FRM1N*20,0 

CALL.  PLTAUtXST  ART,. 5_..QdU2iiiMfc_G4 ,2) 

XSTART«FRM1N*1.0 

CALL  PLTALiXSTARTr,7.ai.ib25jJP0L<KK),5) 

CALL  PLTALCXSTART.Jl, .0.25i INTBBiA) 

CALL  PLTALlXSTART.l?.  ,0.25,  INTEGRA) 

XDECsFRM  I  N*6  -  9 

CALL  PLTAL(XD6C,2i;,(5L25i.UTfeGX,i) 

CALL  PLTAKXOEC.IV,  ,0.25 1  I  NTfcGX  .  1 ) 

XDECuFRMl N*S ,  Q  . . . 

CALL  PLTAL(XDEC,21. ,0.25, !NEAD(KKK),3) 

CALL  PLTAKXDEC,  19;,  0.25,  LFRkQ(Kli.2) 

SELECT  PRINTER  ..MOM  ..  .  . 

CALL  PLTPT<2.,0.,U) _ 

Y0AT.19, 

SYM0«16 . 

PRINT  RUN  INFORMATION  _ 

PRINT  22,FR6Q,NTHD,NRUN,NREC*ALTKF,NDATE,NTIM6,P0HE,NP0L 
22  FORMAT(lHl.//.5X.ltlHFREOUiNCY«,Fl0.4,3HMCS,5X,12HTRUE  MOOING,  15, 
C5X.7HRUN  NO., I3,/,5X,8HRECEIVER,I3,5X,8HALTITUD6,F5.1,5HKFEET,5X. 
C4H0AT6, I8.4X,4HTIMg, I8,4X,6HPOmER«,F10.2.2HDB,7X,3HPOL. M) 

QO  TO  707  . 

7061  PRINT  33.FREQ,NTHQ,NRUN.NREC»ALTKF,NDATb,NTIME,P0WE,NP0L 

33  FORMAT<//,15X.F10,4.20X, J5,12X,I3./,13X, I3,13X,F5.1,14X, I8.8X, 

1  I8.10X.F10.2.12X, J4) 

CHANGE  AND  SET  SYMBOL 

707  SYMB«SYMB-1,0 

CALL  PLTPT ( SYMB , 0 . ,12)  . 

PRINT  DATE  INFORMATION 
XSTART«FRMIN*5.0 

YDAT.YDAT-2.0 
CALL  PLTPT ( 2 . i 0 . . 19) 

CALL  PLTPT<XSTART,y0AT*0, 375,1) 

DATE«N0ATE 

xstart*frmin»9.o 

CALL  PLTNU(XSTART.y04T,DATE<l> 

RESET  SYMBOL 

CALL  PLTPT ( SYMB ,  0  , , 12 ) 

SELECT  LINE  MODE 

CALL  PLTPT ( 1 . , 0  ,  < 19 ) 
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c 

C  SELECT  LONG  LINS  MADS 
CALL  PLTPTI5.  .0. »4) 

PLOT  PATTERN 

CALL  PLTPTCXXC1) »  Yy ( 1 ) ,  8  ) 

CALL  PLTPT  <  XX  C  2 ) . Yv<2)»l6) 

DO  709  IP=3,NP0JNTi 
CALL  PLTPT(XX(  IP) ,VY( IP) ill) 

709  CONTINUE 
CALL  PLTPT ( 2  .  .0..19) 

PLOT  DATA  POINTS 

DO  710  I P*1 j  NPO I N T l 
CALL  PLTPT(XX<  IP)  *  y  Y  ( I P ) # 1 > 

710  CONTINUE 
GO  TO  911 

8  READ  INPUT  TAPE  2 • 777 . NPT, (XX ( J ) ,  y  Y  (  I ) ,  I  *1,  NPT > 
GO  TO  911 
10  GO  TO  1 
END 
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PROGRAM  DATAVE 


PROGRAM  DATAVE 

DIMENSION  GG(100),SDfiOO),RR(10Q)«NA(100),R(9,100)»G(S,lOO) 
DIMENSION  XX(100),YY<100),NF!13> ,NANG(3S),NPL<10> 

DIMENSION  XTW(IOO) , YTW(IOO) 

DIMENSION  N(6) «F(3) # NPO I  NT  C  5  > 

COMMON  R,G,NRdNS,NPQINT,8G,*D,RR,NA,K,KK,KKK 

DATAMNF! I). I *1,13) *9, 11, 13, IS, 16, 17, 19, 21, 22, 29. 26, 29, 30) 

C  STARR  HILL 

DATA<(NANG<l),I«l,7)*397,2,7,36,183,186,193> 

DATA( (NANG(I), I*8,20)*352,12,17,26,31,4i>46,178,196, 

1  917,911,919,920) 

C  PANAMA 

DATA( ( NANG( I), la 21, 29 >*396, 1,6, 11, 16) 

C  THULE 

DATA ((NANG( I), I *26, 30) *179 ,184, 189, 194, 199) 

C  ICELAND 

DATA((NANG(I),I*31.39)*290,299,260,269,270) 

C 

DATA ( ( NPL (I)»l*l,10)*l»2, 3,4,9* ft, 7, 8,9,0) 

READ  SEARCH  PARAMETERS 

1  RE>D  11 i KK , K , KKK 
11  FORMAT (312) 

IF(K)  10,10,2 
10  STOP 

2  NRUNSaO 
DO  3  1*1.9 

3  NPOlNT { I )«0 

SEARCH  IS  STARTED 

913  READ  INPUT  TAPE  2, 766, < N< I ? » I *1 , 7) , (F ( I ) , I ■! , 3 ) 

766  FORMAT (12,11,16,13,12, I«. 12 ,3E16,9) 

NFO«F (2 ) 

IF(NF0.60,NF (K) , AND«N(4 ) .EO, NANG (KKK) , AND  ,N (2) , EO.NPL (KK ) )f ,  4 

4  IF(N(1) ,EQ.99)  6,9 

search  is  ended— ready  for  averaging 

6  CALL  AVERAGE 
REWIND  2 

READY  FOR  NEXT  PARAMETER 
80  TO  1 

9  READ  INPUT  TAPE  2, 777.NTW, (XTW( I I » YTWtf I ) , 1*1 ,NTW) 

777  FORMAT! 12 ,8£16 , 9/ ( 6fil6 , 9 ) ) 

00  TO  913 

8  READ  INPUT  TAPE  2,777,NPT , (XX( 1 ) , YY( I ) • l«l,NPT) 

C 

C  RECORD  DATA  IN  G  AND  R  ARRAYS 
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NRUNS*NRUNS*1 
NPOINT(NRUNS)*NPT 
DO  7  1 .1 >  NPT 
0 ( NRUNS » I ) *XX ( 1 ) 

IMNRUNS* !)iVY( I) 

7  CONTINUE 
C 

C  continue  SEARCH 

C 

00  TO  913 
END 

SUBROUTINE  AVERA9E 

COMMON  R,  G.  NCURVEi  NP0INT1,  GG,  SD,  RR.  NA.K9,  KK,  KKK* 
DIMENSION  SD(IOO),  RR(IQO),  0A<4) 

DIMENSION  R(S.100)i  Q (9 . 100 ) , NPO I ^T1 ( 5 ) »  NA<100),  GQ(lnO) 
DIMENSION  NF(13},NANG<35jil£QkUJU- 

DATA ( ( NF ( I ) .  l.l,  13) .9, 11, 13',  15.16,17. 19, 21. 22, 25, 2  a,  29, 30) 

C  STARR  HILL 

DATA! (NANQ(I)« I «1. 7) .357,2 .7, 36, 163,188,193) 

DATA ( ( NANG ( I ) , I «*, 20) ■352,12, 17,26, 32,41,4* ,17*. 198, 

1  917,918,919,920) 

C  PANAMA 

DATA({ NANG (I), 1.21,25). 356,1, 6, 11, 16)  , 

C  THULE 

DATA ( (NANG ( I ) , I .26 , 30 ) .179 , 184 , 189 , 194, 199 ) 

C  ICELAND 

DATA ( (NANG ( I )  ,  1  «3l , 35 ) «2>0 , 2*5 , 26_fl_» 265 , 2 70  ) 

DAT  A  C  ( I  POL  <  I ) ,  I»11iO).5HSH.MHi5HSH-HV,5MSH-VH,5HSH-W,  5HPAN-V, 
1  5HPAN-X.5HTHU-V,5HTHU«X,5HlCEcV,5HICE-X) 

TYPE  REAL  LOGIOF 

L0G10F  ( X )  >•  .43429  *  LOGF()() 

C 

C  PUT  CURVES  IN  ORDER-- ( LOWEST  ANGLES  TO  HIGHEST  ANGLES) 

c 

DO  30  N  s  1,  NCURVE 
NP  .  NPO 1  NT 1 ( N ) 

IF ( R ( N, 1 )  -  R(NiNP) )  30,  31,  31 
31  NP2  ■  NP/2 

DO  12  I  >  1,  NP2 
NPP  ■  NP  -  I  *  1 
TEMP  .  R ( N , I ) 

R ( N,  I )  "  R ( N , NPP ) 

R ( N* NPP )  .  TEMP 
TEMP  «  G ( N , I > 

G(N» I )  ■  G(N«  NPP) 

G ( N , NPP )  a  TEMP 
12  CONTINUE 
30  CONTINUE 
C 

C  FOR  ONE  CURVE 
C 

IF(NCURVE.LE.l)  101,103 
101  NNN.NPOlNTl(l) 
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NPOINTaNNN 
00  102  IaliNNN 

|D( I )>0 i 

QQ( 1 )>G(1«  1 ) 

NA( I )al 
RR< l >«R(1 *  I ) 

102  CONTINUE 
00  TO  10 

103  CONTINUE 

FIND  THE  LONGEST  CURVE 
NP0INT  a  0 

00  50  II  «  If  NCURVE 
IF (NP0INTK 1 1 J-NP0INT)  50*  50.  51 
51  NPOJNT  ■  NPOINTIUH 
NC>(  I 

50  CONTINUE 

THIS  SECTION  (THROUGH  STATEMENT  t  IS  REPEATED  FOR  EACH 

DO  *  J  ■  1.  NPOINT 

SET  UP  RANGE  ARRAY 

RR(J>  •  R(NCiJ) 

RRA  ■  RR(J) 

CALCULATE  GAINS  AT  THE  ANGLES 

DO  14  I  «  1*  NCURVE 
NP  a  NPOINTl(l) 

I F ( R C I .NPJ-RRA)  20.  22.  22 
a 2  IF(Rd.l)  -  RRA)  21a  21*  20 

20  6 A ( 1 >  a  100. 

GO  TO  14 

21  DO  1  K  a  1.  NP 

IF(R{ I *K)-RRA )  1,  2,  9 

2  GA{ I )  a  0(1, K) 

00  TO  14 

3  K1  •  K-l 

0  A ( I )  a  0  ( I ,  K )  -  (G(I*K)«Q(I*K1))*(R(I *K)«RRA)/(R( I *K) 
00  TO  14 
1  CONTINUE 
14  CONTINUE 

CALCULATE  AVERAGE  0A1N 

GG< J5  a  0,  3  S  a  0.  S  NA(J)  •  0 

DO  5  L«i. NCURVE 
IF(GA(L)-50, )  7,  7,  5 
7  NA(J)  a  nA(J)  *  1 


POINT 


R( I iKl) ) 
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POWER  ■  (5  A  C  L  )  / 1 0  » 

G  A  ( L  >  c  10, ••POWER 
GG(J)  ■  GG ( J )  *  0 A ( L ) 

5  CONTINUE 
GO  (  J )  «  GG(  J  ),":*(  J) 

GGREGzGG ( J ) 

GQ(  j)  *  10.  •  LOG10F  (  G  tt  <  J  > ) 

calculate  standard  deviation 

DO  8  L  •  1>  NCURV6 
I F ( GA i L  >  -  SO,)  40  ,  40,  I 
40  8  ■  S  ♦  { G A  <  L  > "  GGREG ) **2 
9  CONTINUE 
S  »  S/NA<J) 

S  ■  SQRTF(S) 

60 ( J)»(4,3429«S)/QGR§G 

READY  FOR  NEXT  POINT 

6  CONTINUE 
PRINT  INFORMATION 

19  CONTINUE 

PRINT  98  »  NF ( K9 ) «  IPOt-<KK>  »NaNI»(KKK) 

99  f  ORMAT<//,5X»10HrSBJ3U6NCY?,IRl5XilBHSl.Tfe»POLARI2ATUUI»#2X,A5.5X» 

1  6WANgLE«,I4) 

PRINT  999 

999  FORMA T(12X,8HAVE  G A  I N , 15X , 5H ANGLE , 12X , 8HSTD . DE V . , 6X . 6HP0 1  NTS ) 
PRINT  99 , ( GG ( I ) , RR ( I ) , SD ( I ) , N A ( I ) , I »1 , NPO I  NT ) 

99  FORMAT ( 3F20 . 4 ,  10x.  J2> 

READY  TO  PLOT  INFORMATION 

908  CALL  PL0T(GG,RR,K9,KK,KKK,NP0INT) 

108  RETURN 
END 

SUBROUTINE  PLOT < XX  ,  YY , K , KK , KKK , NPO I  NT! ) 

DIMENSION  XX(IOO) , YY  C 100  » 

DIMENSION  I  FREQ (13) i I  HE AD { 35 )  ,  I  POL ( 1 0 > 

DATA((1FREQ(I),Ir1i13)b2H9  . 2MU , 2H13 , 2H15 , 2H18 , 2H17 , 2H19, 2m21 , 

C  2H22»2H25#2H26»2H29,2H30) 

C  STARR  HILL 

DATA( ( 1  HE AD ( I ),  | «1 , 7 ) »3H357 , 3H0Q2 , 3H0Q7 , 3H03* , 3M183 , 3H188 , 3M1P3 ) 
DATA ( (  IHBADd),  I»8,20)»3M352.3HQl2.3H017l3MQ2613H03l,3H04li3H(i**, 
1  3k178,3h198,3H917,3H918,3m919,3W920) 

C  PAN,1  MA 

DATA ( (  I  HE AD ( I ) , I «21 , 25 ) ■3H3S6 , 3W0Q1 , 3M006 , 3M011 , 3H0J6) 

C  THULE 

DATA ( ( 1  HEAD ( I ) , I«26,30»«3H179,3H184,3H189,3H194,3H199) 

C  ICELAND 


125 


DATA <UHEAD(  I  >.1*31,39 )«3H290,3H255,3M2*0,3H249,3h270) 

C 

DATA( ( I  POL ( I ) . I*1»10)*5H8H-MM,SHSM-HV,5WSH-VM,SMSH-VV,9HPAN.V, 
C  SHPAN-X.9HTHU*V  »9HTHU"X, 3HI  *E-v,!>HICi-x) 

DATA  (NITART-0) 

IF  INSTART  .EQ.  0)  <04 • *02 
404  CONTlNUl 
NITART«1 
NFREQP«0 

CALL  PLTPT ( 300  • » 300 , , 2 ) 

CALL  PLTPT(4, ,3, ,13) 

CALL  PLTPT<5. 0,0,0,14) 

•02  IF (RFREOP  .EQ,  K)  604,421 
421  CONTINUE 

CALL  PLTPT ( 0 . • 0 , » 5 ) 

NFREOPaK 

RHIN*XX(1) 

NRMJNaRMIN-9,0 
NXOF8T ■-< ( NRM I N-l ) /4 ) 

XOFST  ■  NXOFST 

CALL  PLTPT(X0FST.1.,3) 

XM0F3T  ■  -X0F8T 

frmin«xmofstm,o 

CALL  PLTSCtFRM INi 4  i,10. ,1,-3. ,3, #8, .1) 

CALL  PLTPT(1. ,0.,10) 

CALL  PLTPT tl • # 0 , ,4) 

INTEQMHHEAQ 
INTE02.IHFREQ 
INTEQ3  ■  6HQEQ 
INTEQ4  ■  SHOa 
INTEQXbIH- 
XSTART«F»HIN-9.0 

call  PLTAL(XST*RT,12i ,0,25, JNTEQ3«3) 

XSTART  •  FRHIN  *  20.0 

CALL  PLTAL(XSTART.-9.0,0,29iINTEQ4,2) 

xstart«frmin*i , 0 

CALL  PLT4L(XSTART,-7,0,0,29, IP0L(KK).9> 

CALL  PLTAL<XSTART.21iQi0.?’iINTE02'4) 
call  PLTALIXSTART ,  19, 0 ,0, ';S, INTEB , 4 ) 

XDEC»FRMIN*6.0 

CALL  PLTALI XDEC,21 , ,0 ,29, INTEQXt 1 ) 

XDEC*FRHIN*8. 0 

CALL  PLTAL(XDEC,2l..0.29,IFREQ(K)  .2) 

HR03«21.0 

SYH0.16, 

•04  HROS«HP08-2,0 
707  SYHS  8  8YMB-1. 

CALL  PLTPTC2.,0.,19) 

CALL  PLTPT(SYM8iO,,12) 

XDEC*FRHIN*9.0 

CALL  PLTPT(XDEC,HPOS*0. 379,1) 

CALL  PLTPTC1. ,0.»19) 
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CALL  PLTPT ( 1 . , 0 , ,4) 

XDEC*FRMIN*6 . 0 

CAL!  PLTAL(XDEC.HPOSi 0.23, 1NTEGX.1) 
XDEl'»FRM  I  N*8 , 0 

CALL  PLTAL(XDEC.HP0S*0.23,  IMEADU'XK)  ,3) 
CALI  PLTPTC5. .0. .4) 

CALI.  PLTPT(XX(1)#YY(1).8) 

CALI,  PLTPT(XX(2)  ,  YY  ( 2 )  ,  16 ) 

DO  'C9  I Pa3 , Npo J NT1 

CALI.  >LTPTCXX(  IP)  «YYCIP).ll) 

709  CONTINUE 

CALL  PLTPT  <  2 . iO, .19) 

CALL  PLTPTISYMB.O, ,12) 

DO  710  IP«i  iNPOINTl 

call  pltptixxi jp>  ,  yyiip)  ,n 

710  CONTINUE 
RETURN 
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